x_av=mean(gaussdata);
s=0;
for i=1:1:size(gaussdata,2)
    s=s+(gaussdata(i)-x_av)*(gaussdata(i)-x_av);
end
s=s/(size(gaussdata,2)-1)